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ABSTRACT 

The nearby transiting planet HD 189733b was observed during three transits with 
the ACS camera of the Hubble Space Telescope in spectroscopic mode. The resulting 
time series of 675 spectra covers the 550-1050 nm range, with a resolution element 
of ~ 8 nm, at extremely high accuracy (signal-to-noise ratio up to 10,000 in 50 nm 
intervals in each individual spectrum). Using these data, we disentangle the effects of 
limb darkening, measurement systematics, and spots on the surface of the host star, to 
calculate the wavelength dependence of the effective transit radius to an accuracy of 
~ 50 km. This constitutes the "transmission spectrum" of the planetary atmosphere. 
It indicates at each wavelength at what height the planetary atmosphere becomes 
opaque to the grazing stellar light during the transit. In this wavelength range, strong 
features due to sodium, potassium and water are predicted by atmosphere models for 
a planet like HD 189733b, but they can be hidden by broad absorption from clouds 
or hazes higher up in the atmosphere. 

We observed an almost featureless transmission spectrum between 550 and 1050 
nm, with no indication of the expected sodium or potassium atomic absorption fea- 
tures. Comparison of our results with the transit radius observed in the near and 
mid- infrared (2-8 /an), and the slope of the spectrum, suggest the presence of a haze 
of sub-micron particles in the upper atmosphere of the planet. 

Key words: planetary systems - methods: data analysis - techniques: spectroscopic 



1 INTRODUCTION 

Transiting planets offer a unique opportunity to study the 
atmospheres of extrasolar gas giants. During the transit, the 
height at which the planetary atmosphere becomes opaque 
to the grazing star light varies with wavelength. Thus the 
transmission spectrum of the planet's atmosphere along its 
limb can be recovered from the wavelength dependence of 
the radius of the planet indicated by the depth of the tran- 
sit. The transit depth change is very small - typically less 
than one part per thousand of the stellar flux - but it 
is accessible with space telescopes for planets with bright 
enough host stars. Transmission spectroscopy provides in- 
formation on the composition and physical properties of the 
planetary atmosphere close to the limb. Several theoretical 
studi es have predicted transmitted spectra of hot Jupiters 



studi es nave predicted t ransmitted spe ctra ot not Ju piters 
fe.g. ISeager fc Sasselovl l200d : iBrownl 1200 ll : iBarmanl (20071 : 
iTinetti et al.l 120071 ). In the 550-1050 nm wavelength range 



covered by the present study, the models predict strong 
atomic lines at 589 nm from sodium and 769 nm from potas- 
sium, and several molecular bands due to water. Further 
into the red, the model spectra are dominated by molecu- 
lar bands of water, methane, ammonia and carbon monox- 
ide. The largest degree of uncertainty in the models is be- 
lieved to be governed by the possible presence of condensates 
forming hazes or clouds, that can block the stellar flux at 
any height in the atmosphere. Hazes and clouds produce 
featureless continuum absorption in the spectrum, with a 
shape depending on the size distribution of the particles. 
In close-in gas giants (T cq = 1000 — 150 K), silicates and 
iron c ompounds can form solid particles l|Seager fc Sasselovl 
l200d : iFortnevI 120051 ). Condensates can also arise from the 
formation of compounds by non-equilibrium p hotochemistry 
in th e strongly irradiated upper atmosphere (jFortnev et al.l 
2003). Strong signatures of clouds and hazes are observed in 
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the spectr a of brown dwarfs with similar atmospheric tem- 
peratures l|Helling fc Woitkell2006l ). 

All planets of the Solar System with a thick atmo- 
sphere have spectral properties strongly affected by hazes 
and clouds. Venus, Titan and Triton have opaque atmo- 
spheres in the visible, that become transparent to much 
deeper layers in the mid-infrared, where the wavelength ex- 
ceeds the typical size of the condensate particles by a suffi- 
cient factor. 

The atmospheric properties of transiting extrasolar gi- 
ant planets are inc reasingly constrained by new observa- 
tions . HD 209458 (ICharbonneau et all |2000| ; iHenrv et af] 
I2000T I and HD 189733 jBouchv et al.1 12005T ) are the most 



scrutinized systems because their bright host star (V ~ 7.7 
mag) and large radius allow precise enough spectroscopic 
measurements. The first positive application of transmis- 
sion spectroscopy was the detection of the Nal absorption 
signature by the Hubble S pace Telescope (HST) in the at- 
mosphere of HD 209458b (|Charbonneau et al.|[2002T ). This 
signature is much weaker than expect ed, which was inter - 
preted as due to the presence of clouds (|Fortnev et al.| [2003). 
photoionization (iBarmanl 1200711 or condensation of sodium 
on the night side (jlro et al.l 12005). So far all ground-based 



upper limits (e.g. Moutou et al. 
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2005 



2007). The varying spec- 



tral signature of the Earth's atmosphere can easily mask the 
tiny effect of the differential transit radius of the planet. The 
much higher precision of data obtained from space is critical 
in this search. 

HD 18973 3b is a gas giant p lanet orbiting a K dwarf star 
in 2.2 days (|Bouchv et alj|2005l ). Because the primary star is 
small (0.755 R©), this system offers the advantage of a very 
deep transit signal (2.4 % flux drop). Data obtained by the 
Spitzer Space Telescope in the mid-infrared on this system 
allowed the analysis of the primary transit llKnutson et al.l 
120071 ; iTinetti et a"ill2007l; lEhrenreich et all |2007|). the sec- 



ondary eclipse llDeming et al.l 2005 j_ Grillmair et al.l 120071 ; 
iKnutson etal] 120071 ; ISwain et ail l2007t). and the day-to - 
night flux va r iation along the orbit I Knutson et alj 120071 ) . 
ITinetti et all |2007l ) present a complete model of the trans- 
mission spectrum fitted to the mid-infrared data (Spitzer 
bands at 3.6, 5.8 and 8 /im), and predict very strong features 
at shorter wavelengths for an atmosphere without conden- 
sates. An HST lightcurve of HD 189733 in the 550-1050 nm 
range has been obtained over three transits with the Ad- 
vanced Camera for Surve ys (ACS) and analysed in white 
light bv lPont etail (|2007l . hereafter "Paper I"). This data 
set provides the opportunity to test the prediction of Tinetti 
et al. in the 600-1000 nm range, and in particular to explore 
the presence of water absorption and the transparency of 
the atmosphere of the planet. The signal-to-noise ratio of 
the ACS spectra is much higher than the Spizer mid-infrared 
lightcurves because of the larger size of HST and the higher 
stellar flux in the visible. The ACS data is precise enough to 
distinguish the signature of a clear atmosphere with water 
and alkali metal absorption features from that of clouds or 
hazes blocking the stellar flux higher up in the atmosphere. 

In this paper, we present the spectroscopic analysis of 
the ACS/HST data set, in order to constrain the trans- 



mission spectrum of the planetary atmosphere from 550 to 
1050 nm. 



2 THE DATA SET 

The data set analysed in this paper consists of three time 
series of ACS low-resolution spectra covering five HST or- 
bits each, during three transits of the HD 189733 system. A 
total of 675 spectra were measured, covering the 550-1050 
nm wavelength range, with a resolution element of about 8 
nm. The accuracy on the relative intensity integrated over 
all wavelengths is 7 x 10 -5 on individual exposures, and 

3 x 10 -5 on 10-minute averages, including correlated non- 
photon sources of noise. The data is presented in detail in 
Paper I, together with the analysis of the lightcurve ob- 
tained by integrating the first-order ACS spectra over all 
wavelengths, and decorrelating for instrumental systematics 
related to changes in the telescope and detector conditions. 
Figure [T] shows this "white" lightcurve (i.e. integrated over 
all wavelengths) of our ACS/HST data. The only significant 
departures from a transit curve model are observed near the 
middle of the transit, and attributed to the occultation of 
starspots by the planet. In the present paper, we analyse 
the chromatic information in the spectroscopic time series, 
with the objective of measuring the transmission spectrum 
of the planetary atmosphere during the transits. We use the 
50-nm-wide binning of the spectrum described in Paper I, 
i.e. our primary data set consists of ten time series, each cov- 
ering a 50-nm wavelength interval from 550 nm to 1050 nm. 
These are obtained by restricting the extraction to specific 
column ranges. We also repeated the reduction using the 
finest sampling possible, 1-pixel wide intervals, which corre- 
sponds to about 3-nm bands in wavelength. This produces 
173 passbands from 537 to 1069 nm. The chromatic data 
is decorrelated from instrumental systematics in the same 
manner as the white lightcurve. The fine-grain data set is 
useful to examine the structure of narrow spectral features, 
namely the sodium and potassium atomic lines. The coarser- 
grain extraction (50-nm bands) is more robust against errors 
in the decorrelation of systematics, and we use it to measure 
the overall shape of the spectrum and the effect of molecular 
bands. 



3 DATA ANALYSIS 
3.1 Overall strategy 

The depth of the dimming due to a planetary transit is di- 
rectly proportional to the projected surface of the planet 
relative to the star, which grows with the square of the 
planetary radius. In a monochromatic lightcurve, the transit 
depth measures the square of the radius at which the plan- 
etary atmosphere becomes opaque to the grazing star light 
for a specific wavelength. The transmission spectrum of the 
planet is therefore obtained by fitting the same transit curve 
model for different wavelengths, leaving the planetary radius 
as a free parameter. In wavelengths where the absorption of 
stellar photons passing through the planetary atmosphere 
is more efficient, the planet appears to have a larger radius. 
The principle is simple enough, but requires data of exquisite 
quality, because the signal represents only a tiny fraction of 
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Figure 1. The HST/ACS lightcurve integrated over all wave- 
lengths. In the middle of the first transit, the planet occults a 
large starspot, producing a steep rise in intensity at the end of the 
second HST orbit. Smaller structures on the star are also crossed 
during the second transit. See Paper 1 for a detailed analysis. 



the signal of the transit (typically 10~ 2 to 10~ 3 ), which is 
itself a small modulation of the star's lightcurve. 

Among known transiting planetary systems, HD 189733 
is the system where the transmission spectrum is expected 
to be most accessible to observations among known transit- 
ing planets, because the radius ratio is large, the planet has 
a hot H2 atmosphere with a large scale height, and the star is 
bright. Signatures of radius change in the transmission spec- 
trum are expected to amount to 1CP 3 — 1CP 4 of the stellar 
intensity. However, the case of HD 189733 is complicated 
by the fact that the host star is variable. It is active at the 
few percent level in visible photometry, with obvious mod- 
ulations of the lightcurve due to rapidly- varying starspots , 
rotating with the st ar in about 12 days (|Croll et alj 120071 ; 
iHenrv fc Wirmll2007l ). The presence of spots has two effects 
on the measurement of the transmission spectrum: (1) the 
occultation of spots by the planet during the transits in- 
troduces features in the transit lightcurve (2) unocculted 
spots on the surface of the star modify the relation be- 
tween the transit depth and the radius ratio. Both effects 
are wavelength-dependent, and must therefore be taken into 
account not only as a global correction to the transmission 
spectrum, but as a passband-to-passband correction. 

The challenge is to disentangle properly the four 



following causes for the variability of the spectroscopic time 
series: 

- Main transit parameters : period, epoch, system scale 
(a/R), orbital inclination angle 

- Telescope/detector systematics 

- Limb darkening 

- Spots and variability 

-Wavelength dependence of the radius ratio (transmission 
spectrum) 

Thanks to the extremely high signal-to-noise ratio of the 
ACS data, a sufficient precision on the transmission spec- 
trum can be achieved to detect the main spectral features 
predicted by planetary atmosphere models. The signal-to- 
noise ratio on individual data points in each 50-nm passband 
is in the 4000:1 - 10000:1 range. With 675 images and 10 
passbands, there is an enormous amount of information in 
the data itself. 

In the following four subsections we consider the causes 
of variability one by one, first explaining how we tackle them, 
and then estimating what level of systematic uncertainties 
they are likely to leave in the measurement of the transmis- 
sion spectrum of HD 189733b. 

3.2 Main transit parameters 

The orbital period, phase and inclination angle, and the 
system scale, do not depend on the wavelength. We fix 
the values of these parameters to those derived in Paper I 
from the lightcurve integrated over all wavelengths ( "white" 
lightcurve) . 

The effect of the uncertainties on these parameters can 
be evaluated by moving them by their ±1-<t error interval 
and repeating the whole analysis with these new values. The 
result of this procedure is shown on the upper left panel of 
Fig. [5] Changing the transit epoch and period has little ef- 
fect on the transmission spectrum. Modifying the inclination 
angle has more impact, both on the absolute value and on 
the slope, because the run of limb darkening along the chord 
of the projected path of the planet on the star changes. In all 
cases, the uncertainties on the main transit parameters pro- 
duce changes with a very smooth wavelength dependence, 
so that they cannot introduce or mask narrow features in 
the transmission spectrum. 

3.3 Telescope and detector systematics 

Systematic flux variations in the ACS data were removed 
with decorrelation using six vectors: the x and y position 
on the chip, the spectrum width, angle, the HST orbital 
phase, and linear time. The details are given in Paper I. 
The passband-integrated intensities were decorrelated in the 
same way as the white lightcurve, on each passband inde- 
pendently. We did not attempt to take into account the co- 
variance between the corrections of neighbouring passbands. 

Residual systematics can be estimated from the level 
of correlated noise in the extensive part of the time series 
that was taken outside the transit. The top right panel of 
Fig-[l]shows the amplitude of the residual system atics in the 
ten pa ssbands that we use, using the method of lPont et al.l 
(2006) to estimate the effect of red noise on the planetary 



4 F. Pont et al 



— 400 
B 

u 300 

tic 

jS 200 

o 

.2 ioo 

tO 

£ o 
S-ioo 



- 1 1 1 1 1 1 1 


111111111111111 


1 1 - 


! t 


Main 


1 

p ararnsLers 




6 


- 






m 


- 






i — i 


— 






CD 
















G 








to 
















o 








OT 
















CO 




□ □ □ 


3 : 


Sh 


XXX 


X X x x X X 




















1/5 








e 








to 










"1 1 1 1 1 1 1 


II 


i - 


E-i 



600 700 800 900 1000 
X[nm] 



600 700 800 900 1000 
X[nm] 



400 
300 



100 


100 



i i i i I i i i i I i i i i I i i i i I i 
Telescope/ detector 




600 700 800 900 1000 
A[nm] 




i i i i I i i i i I i i i i I i i i i I i 
Unocculted spots 



□ □ 



Pi i i I 



600 700 800 900 1000 
A [ n m ] 



Figure 2. Changes in the transmission spectrum results caused by the uncertainties in the different steps of the analysis. The vertical 
scale shows the change in planet transit radius. The dotted bins in all panels show the photon shot noise. Top left: effect of moving the 
transit timing (crosses), radius ratio (open squares) and inclination angle (solid squares) by the 1-cr uncertainty interval in P aper 1. Top 
right : effect of residual red noise left over from the decorrelation of instrument /detector parameters, using the method of lPont et al.l 
(2006). Bottom left: effect of shifting the limb darkening coefficients by one passband. Bottom right: effect of adding unocculted 
spots absorbing 1% of the stellar flux (open squares), and of changing the mean temperature of the spots from 4000 K to 3500 K (solid 
squares) . 



parameters. They have an amplitude comparable to the pho- 
ton noise, except at the edges. The shortest and longest pass- 
bands are more sensitive to systematics because the actual 
flux on the detector drops sharply on the edges of the spec- 
trum, which amplifies the effect of telescope pointing drifts 
and focus changes. 

Three different transits are covered by our data, so that 
the instrumental systematics can also be estimated from the 
differences between the transmission spectrum measured in 
the three visits. The systematics are not correlated between 
the visits, because they were taken at three very different 



positions in terms of HST configuration. Of course, this 
estimate will be an upper limit, because other factors af- 
fect the comparison between the three transits (namely the 
starspots, see below), and the second and third visits have 
fewer points during the transit. The comparison between the 
results from the three visits is discussed in Section [4] 



3.4 Limb darkening 

Limb darkening has an important effect on the transit curve 
shape in the wavelength range covered by the data, and 
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Passband [nm] 


cl 


c2 


c3 


c4 


550 


- 600 


0.4621 


-0.2003 


0.9450 


-0.4045 


600 


- 650 


0.5148 


-0.2774 


0.9429 


-0.4033 


650 


- 700 


0.5620 


-0.3531 


0.9524 


-0.4096 


700 


- 750 


0.5932 


-0.4183 


0.9488 


-0.4003 


750 


- 800 


0.6231 


-0.4808 


0.9672 


-0.4078 


800 


- 850 


0.6390 


-0.5260 


0.9695 


-0.4039 


850 


- 900 


0.6529 


-0.5366 


0.9290 


-0.3886 


900 


- 950 


0.6538 


-0.5508 


0.9273 


-0.3867 


950 


- 1000 


0.6475 


-0.5439 


0.9104 


-0.3837 


1000 


- 1050 


0.6303 


-0.4967 


0.8538 


-0.3696 



Table 1. Flux- weighted limb darkening parameters for our pass- 
bands calculated from a Kurucz model atmosphere. 



is strongly coupled with the transit depth. Fortunately, 
limb-darkening is a well-understood effect and can be mod- 
elled with high accuracy for a K dwarf in this range of 
wavelengths. We used a synthetic stellar spectrum calcu- 
lated by R. Kurucz for T eff = 5050 K logg = 4.53 and 
[Fe/H]=-0.03 (from iBouchv et aL|[2005l rl. to compute the 
intensity in small wavelength intervals and at 17 radial posi- 
tions spanning the region from the center out to the edge of 
the star. We then calculated the flux-weighted average in- 
tensity distribution across the surface of the star in each 
of our passbands, and fit this intensity distribution with 
the four-para meter non-linear limb-darkening law given in 
IClaretl (|200(t ) (J(/x) = 1 - £* =1 c ?l (l - M n/2 )) to determine 
the limb-darkening coefficients for each passband. The cor- 
responding limb-darkening coefficients are listed in Table [T] 
To estimate the potential systematic uncertainties in- 
troduced on the final spectra by the choice of limb- darkening 
coefficients, we computed the spectrum using for each pass- 
band the limb-darkening coefficients appropriate for the next 
passband. The bottom left panel of Fig. [2] shows the dif- 
ference between the resulting spectrum and the one using 
the limb-darkening parameters of Table [1] The effect of this 
rather large simulated uncertainty on the limb-darkening 
parameters is to shift the whole spectrum by an amount 
comparable with the photon noise. As expected, a change 
in limb-darkening parameters affects the overall slope of the 
spectrum, via the covariance between limb darkening and 
transit depth. Spectral features, however, are not affected, 
because the intensity and shape of limb darkening varies 
rather smoothly in the A = 600 — 1000 nm range at that 
spectral resolution. 

3.5 Starspots 

The presence of starspots on the surface of HD 189733 affects 
the measurement of the planetary transmission spectrum in 
two ways: 

(a) - spots occulted by the planet modify the shape of the 
transit lightcurve. The modification is chromatic, since the 
occulted spots have a different spectrum from the rest of the 
star 

(b) - spots not occulted by the planet, but visible to the ob- 
server, diminish the brightness of the unocculted part of the 

1 The synthetic spectrum for HD 189733 is available on the web- 
site of R. Kurucz, |http://kurucz. haxv ard.edu/stars/HD189733 



star compared to a spot-free surface, and therefore slightly 
increase the depth of the transit. This is also a chromatic 
effect. 

One big starspot and a series of small spots are occulted 
during our time series, and extensive ground-based and 
space-based follow-up constrains the distribution of spots on 
the whole surface during th e HST observations (|Croll et al.l 
120071 ; iHenrv fc Winnll2007l . see Paper I). There is therefore 
enough information available to correct for the effect of the 
spots. 

(a) Spots occulted by the planet 

The planet occulting a spot on the surface of the star 
during the transit produces a rise in the lightcurve, because 
the spot is cooler and fainter than the rest of the star. Hiding 
the spot therefore increases the overall luminosity of the part 
of the star not hidden by the planet. This effect varies with 
wavelength, according to the difference in surface brightness 
between the spot and the rest of the stellar surface, i.e. the 
effective temperature of the spot. To first order, the effect 
behaves like the difference of the spectra of two stellar pho- 
tospheres at different temperatures. 

In our reference solution, we simply do not use the part 
of the time series visibly affected by starspots. This in fact 
includes more that half of the in-transit data (see Paper I), 
but there is a sufficient stretch of data in the first HST visit 
without apparent spot effects to cover all phases of the first 
half of the transit. Since the transit model is symmetrical 
around the transit mid-time, this sequence is sufficient to 
constrain the chromatic dependence of the transit depth (i.e. 
the transmission spectrum of the planet). 

In order to check the robustness of the reference solu- 
tion, we use the parts of the time series affected by spots 
to compute independent estimates of the transmission spec- 
trum. To first order, the occultation of the large spot - still 
ongoing at the end of the HST sequence - corresponds to 
a linear flux rise. We therefore fit a transit model with a 
chromatic linear rise to the time series. The free parameters 
are the transit depth in each passband, the time of contact 
with the spot, and the slope of the rise in each passband. 
The results in terms of slope vs. passband for the spot are 
shown by the circles in Fig. [3] 

This procedure has an important side benefit: it pro- 
vides the spectral dependence of the effect of the occulta- 
tion of the spot. We can compare the results with the effect 
expected from occulting a stellar surface of lower tempera- 
ture against the background of a T = 5000 K stellar pho- 
tosphere. The expectations for synthetic Kurucz stellar at- 
mosphere models are shown in Fig. [3] for spot temperatures 
of 4500 K, 4000 K and 3500 K, assuming T cft = 5000 K for 
the star as a whole. The intensity rise slope is not a smooth 
function of wavelength, because of the different strength of 
temperature-dependent absorption features on the photo- 
sphere and in the spot. From Fig. [3] we can estimate the ef- 
fective temperature within the spot. The observed chromatic 
effect of the spot corresponds to the expectations from the 
occultation of a cooler surface with effective temperature 
between 4500 K and 4000 K . A lower spot temperature 
of 3500 K predicts a stronger difference than observed in 
the blue, due to the quasi-disapearance of the blue part of 
the spectrum under molecular bands at these temperatures. 
This range of temperature is consistent with what is known 
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Figure 3. Slope of the relative intensity rise caused by the planet 
occulting a starspot during the first transit covered by our data, 
as a function of the central wavelength of the 50-nm spectral bins. 
The lines show the signals expected for spots with temperatures 
of 4500 K (dashed), 4000 K (solid) and 3500 K (dotted), with 
square symbols marking the integration over our passbands, for 
each model. 



for spots on solar-type stars: large starspots ar e 500-1000 K 
coole r than the photosphere in active stars dFrasca et al.l 
l2005h . 

The effect of smaller starspots occulted during the sec- 
ond HST visit is too small to be modelled in this way. We 
can use the model of the large spot, however, to correct the 
effect of the small spots differentially. We assume that the 
small spots have a similar temperature as the large one. We 
compute, on the "white" lightcurve, the difference between 
the time series affected by spots, and the theoretical tran- 
sit curve with the parameters of Paper I. We then remove 
this difference from the coloured lightcurve, multiplied by 
the passband dependence obtained for the large starspot. 
Namely, if fa is the slope of the rise due to the large spot in 
passband i, f the slope in the white lightcurve, the vector 
AF = F — F c is the flux difference between the observed 
white lightcurve and the transit model, then the corrected 
lightcurve for passband i is: 

F[ = F - AF x fi/f 

Obviously, the assumption that the spots have identical tem- 
perature can be questioned. In our case AF is of the order of 
5 ■ 10~ 4 at most, and fi/f varies between 0.8 and 1.2. Con- 
sequently the uncertainty on the exact spot temperature is 
a small effect (<< 2 • 10~ 4 ). 

(b) Spots not occulted by the planet 

Even spots not occulted by the planet affect the transit 
signal. Spots make the unocculted part of the star dimmer 
than a spot-free surface, and therefore increase the depth of 
the transit and the calculated radius ratio. This also depends 
on wavelength, since the dimming due to spots depends on 
the difference of their spectrum with that of the rest of the 



stellar surface. Such effects are very tiny (< 10 -4 on the 
wavelength dependence of the transit depth even for an ac- 
tive star like HD 189733), but they need to be taken into 
account at the level of accuracy reached by our data. 

In Paper I, we accounted for this effect differentially, 
using the variations of the total flux due to spot extinc- 
tion between the three transits covered by the data. The 
flux variations are measured by the absolute calibration of 
the HST data itself, and by the simultaneous ground-based 
monitoring of the star. Figure [4] places the three HST mea- 
surements in the con text of the long-term monitoring by 
I Henry fc Winn! l|2007l ) and the observin g campaign of the 
MOST space telescope (|Croll et al.ll2007l ). The ground-based 
and MOST data are taken at similar wavelengths, centered 
near 560 nm. To account for the much redder position of the 
HST spectra, the intensity scale was expanded by a factor 
2. This factor was found by modelling the spectral distri- 
bution of the effect of spots at 4000 K on a stellar surface 
at 5000 K. In the figure, one complete cycle of the MOST 
data is repeated six times before and after the actual cover- 
age. This extrapolation connects smoothly with the ground- 
based data of the previous and subsequent season both in 
amplitude and phase, which suggests that the pattern of 
spots has remained relatively stable. In contrast, just before 
our first observations, the ground-based photometry indi- 
cates a change of behaviour of the spot pattern in less than 
one rotation cycle. The three absolute flux measurements 
from the HST data, showing variation of +5 x 10~ 3 in flux 
between the first and second visit, and +7 x 10~ 3 between 
the first and third visit (Paper I), are in agreement with this 
picture. 

For the chromatic data, we need not only an estimation 
of the flux difference caused by star spots moving in and out 
of view during the visits, but also of the total flux dimming 
caused by all unocculted spots visible on the stellar surface 
at the time of the transit, because these spots diminish the 
flux of the star throughout the transit, and therefore the 
transit radius calculated from the lightcurve is artificially 
increased. To account for this effect, we have to estimate the 
"background level" of starspots, i.e. the minimum quantity 
of dimming due to spots whatever the viewing angle. 

Croll et al. (in prep.) present a spot model for the 
MOST lightcurve covering two rotation cycles of HD 189733 
a short time after our sequence, and find that a good fit can 
be produced by a two-spot model, with a negligible back- 
ground spot value. The third HST visit is near the epoch 
when the two large spots are almost entirely hidden. How- 
ever, as pointed out by these authors, more complex models 
with numerous, smaller spots, and time variations in spot 
temperature and size, or differential rotation velocities for 
spots as a function of the latitude on the star, could also 
account for the observed lightcurve. In particular, a higher 
angle between the stellar spin and the orbital axis of the 
planet is required in the Croll et al. model than observed by 
IWinn et all (|2007f ) using the Rossiter-Mc Laughlin effect. A 
model with more than two spots would be needed to rec- 
oncile the two pieces of information. Our ACS time series 
tends to indicate that smaller spots are indeed present on 
other parts of the star. 

We have constructed several suites of spot models fitting 
the MOST lightcurve, and have found the following statis- 
tical indications: if a few large spots dominate the signal, 
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Figure 4. Lightcurve of HD189733 over 260 days. Crosses: Ground-based monitoring from iHenrv fc Wirinl d200^ and continuous 
monitoring with MOST ijCroll et al.ll200^ ). Squares: Absolute flux during the three HST visits from Paper I. The zero-points of the three 
lightcurves have been adjusted. The vertical scale is doubled for the HST data, to account for the wavelength dependence of spot-related 
variability. The line shows an extrapolation of the MOST lightcurve over 6 cycles before and after the actual measurements. 



then a background spot level between and 1 % is typical 
(i.e. the spots diminish the intensity by 1% in all viewing 
angles compared to a spot-free surface). For models with a 
large number of spots following a power-law size distribution 
with -1 slope, the typical background spot level is 1-3 %. 

We can also try to estimate the background spot level 
from the HST data itself. Let us consider that the whole 
surface of the star is covered with small and medium spots 
like the ones crossed during the first two HST visits. The oc- 
cultation of these spots causes an average decrease of about 
3 x 10 -4 in the transit lightcurve, or 1.2% of the flux occulted 
by the planet. Extending this to the whole stellar surface di- 
rectly translates into a 1.2% factor for the dimming due to 
background spots if none of the observed photometric vari- 
ation is due to the distribution of these spots, but entirely 
to bigger spots. 

There is a shortcut to control the potential effect of un- 
occulted starspots in this scenario. If we assume that the 
zone of the surface of the star crossed during the second 
visit is typical of the whole stellar surface, then we can fit 
for the effective transit radius ratio without any correction 
for spots. The mean effective radius ratio during this pho- 
tometric sequence will be an unbiased estimator of the real 
radius ratio, since the flux dimming in the part of the star 
occulted by the planet (including the spots) will correspond 
to the dimming by spots of the unocculted part. This will 
be valid in all passbands, provided all the spots have similar 
temperatures. The resulting transmission spectrum will of 
course be noisier, because the spots affect the accuracy of 
the description of the lightcurve by a smooth transit model, 
but the overall level will be correct. The result of using this 
procedure is shown in Fig. [S] (solution 'E'). Note that this 
"worst case scenario" is physically unexpected. Spots in ac- 
tive stars are expected to be large and concentrated at cer- 



tain latitudes, not covering the entire surface l|Granzer et al.l 

l2ood y 

The effect of the uncertainty on the spot temperature 
and on the total amount of unocculted spots is illustrated in 
the bottom right panel of Fig. [21 The black symbols show the 
difference in the transmission spectrum caused by changing 
the mean temperature of the unocculted spots from 4000 
K to 3500 K. The white symbols show the effect of adding 
a constant background of unocculted spots with a spectral 
distribution described by the dots in Fig [3] amounting to a 
dimming of 1% over the whole star. 

Figures [5] and [S] show that the correction for unocculted 
starspots is the dominant source of overall uncertainty on 
the exact shape of the transmission spectrum. It can affect 
its level and global slope at the ~ 400 km level. A change 
in spot temperature has a non-linear effect on the shape of 
the transit spectrum because of the increasing strength of 
molecular bands below 4000 K. 

The consideration of the ACS /HST spectrum in 1-pixel 
bins (~ 3 nm) allows us to place upper limits on the total 
level of background starspots and their temperature. The 
presence of cool starspots on the star introduce spectral fea- 
tures in the apparent radius ratio, because of the temper- 
ature dependence of the sodium line at 589 nm, and the 
presence of molecular absorption bands below 4000 K. The 
absence of any such features in our observed spectra - bar- 
ring a coincidence between the planetary and stellar spectra 
- indicates that the background intensity diminution is not 
higher than about 3 % for spot temperatures near 4000 K, 
and 1 % for temperatures near 3500 K. 

Thus, several independent approaches indicate that the 
background spot level is between zero and 3 %, with lower 
values being more likely, and with typical spot temperatures 
of 4000 K or higher. In our reference solution, we shall as- 
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Figure 5. Planetary transmission spectra obtained from four dif- 
ferent solutions: first orbit without spot (reference solution 'A'; 
solid squares), second orbit with empirical correction of occulted 
spots ('C; diamonds), third orbit ('D'; open squares), and second 
orbit including raw occulted spots signal ('E'; circles). Solution 
'B' is nearly identical to solution A and is not shown. 

sume a value of 1 % for the unocculted spot level during 
visit 1. 



3.6 Wavelength dependence of the radius ratio 

We finally compute the planetary transmission spectrum by 
fitting the ACS time series for the radius ratio, for each 
colour passband. The transit model parameters are fixed 
to those obtained in Paper I from the combined white 
lightcurve, and the limb d arkening coefficients t o the val- 
ues of Table Q] We use the iMandel fc Agoll (|2002h theoreti- 
cal transit shapes. The computed uncertainties take residual 
sy stematics noise in to account following the prescriptions 
of IPont et al.l J200rj). The wavelength-dependent correction 
for unocculted spots is applied a posteriori, passband by 
passband. The spot corrections are taken from the empirical 
method of Section |3. 51 with a 1% dimming from unocculted 
spots. 

We performed five different reductions, labelled A to E. 
In the different solutions we change the HST visits used, the 
way the occulted spots are treated, and the method of de- 
termination of the limb darkening. The results are displayed 
on Fig. [5] in terms of height of the effective transit radius 
as a function of wavelength, relative to an arbitrary fiducial 
level. 

In solution 'A', we use the data from the first HST visit, 
outside of the feature caused by crossing the large starspot. 
This series samples an apparently spot-free region, up to 
the center of the star. It is therefore very sensitive to the 
effective radius ratio. 

The four other solutions are of lower accuracy since they 
use shorter or noisier parts of the time series, but they are 
a very important way to check the overall level of the sum 
of all the sources of systematics, since they are partly or 
entirely independent of each-other. 

In solution 'B', we use the data in Visit 1 including 
the occulted spot. The effect of the spot is assumed to be 
a linear flux increase, with a wavelength dependence con- 
structed from the difference between two synthetic stellar 
spectra from the Kurucz library, with T e g = 5000 K for the 



star and T c g — 4000 K for the spot. The limb-darkening 
coefficients are not fixed to the values of Table [TJ but empir- 
ically determined from the best fit to the data of the three 
HST visits, using a quadratic law. 

In solution 'C, we use the data in Visit 2. The small 
fluctuations due to spots are compensated by assuming the 
same wavelength dependence as the large spot signal ob- 
served in Visit 1. The limb-darkening coefficients are taken 
from Table fT] 

In solution 'D', we use the data in Visit 3, with limb- 
darkening coefficients from Table [1] No correction for oc- 
culted spots is applied. 

Finally, in solution 'E' we use the data in Visit 2 without 
spot correction (the shortcut explained in Section [33]), with 
limb-darkening coefficients from Table [1] 

Overall, the difference between the solutions is consis- 
tent with the level of systematics estimated from the out-of- 
transit fluctuations. Solutions B to E are noisier than solu- 
tion A, as expected from the effects of spots (B, C and E) 
or the presence of less in-transit data (D) . Nevertheless, the 
main characteristics of the resulting spectrum are observed 
in all solutions: the absence of strong feature, a negative 
slope, a rise below 600 nm and near 800 nm, and a decrease 
beyond 1000 nm. Note that solution B and C are entirely in- 
dependent (first vs second HST visit, atmosphere models vs 
fitted limb-darkening parameters, empirical vs model spot 
correction). Solution E confirms that the background level 
of spots should not affect the results fundamentally. Even 
if the stellar surface were entirely covered with small spots, 
the effective radius would decrease by no more than 500 km 
and the global shape would be conserved. 

Solution A is chosen as our 'reference solution', because 
the other solutions are more sensitive to the uncertainties on 
the properties of the spots (solution B, C and E), or to the 
values of limb darkening (solution D). Solution A samples 
all the phases of the transit on a spot-free stellar surface. 
We choose not to combine solutions from the three visits, 
since this would introduce an additional dependence on the 
changing configuration of the unocculted spots. 



4 RESULTS: TRANSMISSION SPECTRUM OF 
THE PLANETARY ATMOSPHERE 

The time series in each 50-nm chromatic passband from 550 
nm to 1050 nm are shown in Fig. [6] The transit depth was 
multiplied by 0.993 (first visit) and 0.998 (second visit), to 
account for the effect of unocculted spots, the sequence of 
the first transit affected by the planet occulting a large spot 
is not included and the second transit is corrected for the 
effect of the small spots (see Section [33] on the effect of star 
spots). 

Table [2] gives our 'reference' transmission spectrum for 
HD 189733b, in terms of transit radius ratio (-Rpi/iistar) in 
ten 50-nm passbands. Figure [7] displays the results in terms 
of height in the atmosphere at which the planet becomes 
opaque, relative to a fiducial height of 82000 km (using 
Rsta,r = O.755i?0 from Paper 1). The boxes indicate, verti- 
cally, the photon-noise error bar, and horizontally the width 
of the passbands. The points are not centered on the boxes 
because the intensity-weighted wavelength means are differ- 
ent from the arithmetic means. Dotted boxes show the am- 
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Figure 6. Phased lightcurves in the ten 50-nm chromatic pass- 
band from 550 nm to 1050 nm, shifted vertically for visibility. The 
lightcurves are corrected for the effect of starspots, occulted and 
unocculted (see text). Longer wavelengths towards the bottom of 
the figure. 



Passband [nm] 


Transit radius ratio 


550 


- 600 


0.156903 ± 0.000095 


600 


- 650 


0.156744 ± 0.000065 


650 


- 700 


0.156552 ± 0.000057 


700 


- 750 


0.156388 ± 0.000059 


750 


- 800 


0.156501 ± 0.000064 


800 


- 850 


0.156210 ± 0.000073 


850 


- 900 


0.156147 ± 0.000081 


900 


- 950 


0.156120 ± 0.000092 


950 - 


- 1000 


0.156097 ± 0.000125 


1000 


- 1050 


0.155716 ± 0.000218 



Table 2. Transmission spectrum of HD 189733b — transit radius 
ratio in ten wavelength passbands, from our reference solution 
with correction for unocculted spots, with photon-noise errors. 
Systematic uncertainties are discussed in the text. 



plitude of the telescope/detector systematics, when larger 
than the photon noise. Note that the scale height of the at- 
mopshere of HD189733b is of the order of 150 km. The un- 
certainties on the transit radius in 50 nm bands correspond 
to ~ 50 km near the central wavelengths. 

The cloud-fre e planetary atmosphere model from 
iTinetti et alli|2007t ). normalized on the Spizer infrared data, 
is plotted for comparison. Strong absorption features from 
sodium (589 nm Na I doublet), potassium (769 nm K I dou- 
blet) and water are predicted. The horizontal dotted line 
shows the transit radius value of iKnutson et al" | (|2007h at 8 
microns, which also corresponds to the highest level of the 
NICMOS transit radius measurements near 2 microns (M. 
Swain, private communication). 

The main overall property of the 0.55 - 1.05 micron 
transmission spectrum of the planet HD 189733 from our 
data is the lack of large features, such as would be expected 




800 1000 
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Figure 7. Atmospheric transmission spectrum of HD 189733b 
in the 600—1000 nm range. The vertical axis shows the height 
in the atmosphere of the effective transit radius of the planet, 
relative to an arbitrary reference level. The points are located at 
the intensity-averaged center of the 50-nm bins in wavelength, for 
our reference solution. The boxes show the extent of the wave- 
length bins horizontally and the photon noise vertically. Dotted 
boxes correspond to the level of instrument systematics when 
larger than the photon noise. T he solid line shows the synthetic 
spectrum of lTinetti et al.l {2007), normalised on the mid-infrared 
Spitzer data, with dots marking averages over 50-nm bins. The 
dotted line shows the infrared radius discussed in the text. 



from sodium, potassium and water in a transparent hot Eb 
atmosphere. The spectrum is relatively flat at the average 
level of Rpx/R* = 0.1564. 

Two more detailed observed properties are noted: a 
smooth slope of ~1 % in transit radius from 600 to 1000 
nm, and a bump in the 750-800 nm passband. These two 
features can be affected by systematics. Nevertheless, they 
are robust in the face of our tests on these systematics. The 
value of the slope, in kilometers of transit radius, is 472 ± 
42 (random) ± 106 (syst) kilometers over a 400 nm inter- 
val. Flattening the slope of the spectrum could be done by 
invoking unocculted spots, but this would require spot lev- 
els that would leave signatures that are not observed in the 
variability pattern and transmission spectrum. 

The mean transit radius in the 600-1000 nm range is 
more than 2000 km higher tha n predicted by the cloud-free 
models of lTinetti et al.1 l|2007l ). and about 1000 km higher 
than the transit radius in the mid infrared. The uncertainty 
on the mean level is ± 16 (random) ± 61 (syst) km. The 
systematic errors are estimated by the quadratic addition of 
all sources of systematics studied in Section [3] 

Slope changes at the 1-a level are observed at both edges 
of the spectra, and correspond to the expected position of 
the sodium line on the blue side, and the edge of a water 
band on the red side. Detailed examination of the pixel-by- 
pixel spectrum shows no indication that these features are 
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due to narrow absorption features. They may be caused by 
the stronger systematics present at both edges. 

No significant signatures of the sodium and potassium 
lines are detected. To quantify the absence of these lines, 
we repeated the extraction in 15-nm bands centered on the 
sodium line at 589 nm and the potassium line at 769 nm, 
compared to a baseline filter of 75 nm around the line posi- 
tions. We measure a flux excess of (—6.9 ±7.9) x 1CP 5 in the 
sodium band (extra absorption during transit would yield a 
positive value here), and (1.6 ±5.2) x 10~ 5 in the potassium 
band. The quoted errors take only the photon noise into ac- 
count. Systematic uncertainties estimated from the actual 
pixel-to-pixel scatter of the spectrum gives ±11.9 x 1CP 5 
and ±21.9 x 10~ 5 for the sodium and potassium bands re- 
spectively. N ote that the 4-sigma det ection of sodium for 
HD 209458 jCharbonneau et all [2000) was with a narrow 
1.2 nm band, their result with a 10 nm band was an in- 
significant (3.1 ±3.6) x 10 -5 dimming during transit. In the 
case of HD 189733, additional unocculted starspots would 
increase the effect of the sodium line in the apparent trans- 
mission spectrum, so that the observed absence of feature is 
an upper limit. The narrow-passband data show no indica- 
tion that the bump in the wide-band spectrum is located at 
the position of the potassium line. However, the wide-band 
feature is visible independently in the spectra from all three 
visits in the HST data. It could therefore be real. 

The data clearly exclude lines wi th the strength pre - 
dicted by the condensate- free models of lTinetti et al.l ((2003). 
This, together with the lower transit radius in the infrared 
and the redward slope in the spectrum, strongly suggests ab- 
sorption of the grazing stellar flux by a haze of condensates 
high in the atmosphere of HD 189733b. 



5 DISCUSSION 

Dust and condensates str ongly affect the emission spec- 
trum of brown dwarfs (e.g. lAllard et al.ll200"ll ; iMarlev et al.l 
2002) with atmospheres similar in composition and tem- 
perature to close-in gas giant planets. Condensates of iron, 
silicates and AI2O3 are thought to b e abundant in brown 
dwarf and hot Jupiter atmospheres (|Fortnevl 120051 ) . With 
the slant viewing geometry of transit transmission spec- 
troscop y, less abundan t compounds could also become im- 
portant. [^rtne3 (|2005h mentions Cr, MnS, Na 2 , Zn, C, KC1 
and NH4H2PO4 with condensation temperatures in the rel- 
evant range. Grains for med by such conden sates are usu- 
ally sub-micron in size l|Helling et alJ kooe). producing a 
continuum absorption peaking in the visible and decreas- 
ing steeply in the near infrared. In contrast to brown dwarf 
atmospheres, hot planet atmospheres are heated from above 
by the irradiation of their host star, so that non -equilibrium 
photo chemical products could also play a role (I Liang et al.l 
2004), like the ozone layer on Earth. 

The observation of an almost featureless spectrum be- 
tween 550 and 1050 nm for the planet HD 189733b is sugges- 
tive of absorption by condensates high in the atmosphere. 
The two strong alkali atomic absorption features expected 
in hot Jupiter atmospheres are either absent or much weaker 
than predicted, which can also be explained by absorption of 



the grazing stellar light by clouds or hazefl before the atomic 
lines can become efficient. The same explanation has been 
proposed to e xplain the weak strength of the sodium line 
in HP 209458 (charbonneau et al.ll200(il ; ISeager fc Sasselovl 
l2000l y 

The ACS, NICMOS and Spitzer data indicate that the 
transit radius in the near infrared is about 1000 km lower 
than in the visible. If confirmed, this difference would be con- 
sistent with the presence of condensates. Particles efficiently 
scatter radiation with wavelengths comparable to their ra- 
dius, then the cross-section drops rapidly. The redward slope 
of the 550-1050 nm transmission spectrum suggests that the 
particles in the haze have sub-micron scales, so that their 
cross-section will be much lower in the near infrared. At- 
mospheres opaque at visible wavelengths but much more 
transparent in the near infrared are known in the Solar Sys- 
tem. On Venus and Titan for instance, condensates block the 
visible light very high in the atmosphere, but much deeper 
layers, or even the planetary surface itself, can be seen in 
some spectral windows beyond 2 microns. Industrial hazes 
on Earth also share similar properties. 

When relating the present results to the reflection and 
emission spectrum of HD 189733b measured during the 
secondary eclipse, several caveats are necessary. First, the 
transmission spectrum is formed solely along the limb of 
the planet, so that the conditions probed are only those of 
the day-night terminator. Second, the convection on close-in 
gas giant planets is expected to operate on a global scale, 
because the synchronous rotation lowers the Coriolis forces 
compared to, for instance, Jupiter. It is therefore concievable 
that the transmission spectrum varies from transit to tran- 
sit, according to "weather" in the upper atmosphere. And 
third, a low opacity sufficient to scatter grazing starlight may 
not represent a dominant contribution with perpendicularly 
incident light. 

With these caveats, the presence of condensates in the 
upper atmosphere could influence the reflection/emission 
spectrum observed during secondary eclipse, through the 
scattering of light as well as by changing the pressure- 
temperature profile. Since most of the flux of the host star 
is emitted in the red, absorption and scattering by the 
haze is likely to modify the temperature profile of the at- 
mosphere by injecting heat at high altitudes. In a similar 
context, a temperature inversion has recently been postu- 
lat ed in the atmos p here of the transiting planet HD 209458b 
by iBurrows et al.l (|2007l) to explain t he em ission spectrum 
recently observed by Knutson et alj (|2007l ) during the sec- 
ondary eclipse. 

As a concluding remark, we note that our results indi- 
cate that the sunset over HD 189733b is red. 
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2 In this context we make no fundamental distinction between 
clouds and haze, which differ only by their opacity. 
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